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INTRODUCTION 


The  subject  of  ore-  and  postprocessing  of  finite  element  data  has 
been  gaining  importance  over  the  last  several  years.  It  is  perhaps 
accurate  to  state  that  research  in  this  area  has'  been  as  active  as 
research  in  the  numerical  method  field.  In  this  context  one  has  to 
distinguish  between  several  functions  associated  with  what  one  may 
call  "auxiliary  finite  element  software".  Pre-processors  perform  such 
functions  as  geometry  definition,  geometry  discretization  (mesh 
generation),  model  display  and  verification,  model  editing,  node  and 
element,  renumber ing  as  well  as  load  and  boundary  condition  definition. 
Operations  such  as  load  and  mass  discretization  are  best  considered  as 
part  of  the  analysis.  Postprocessing  may  be  concerned  with  'result 
interpretation.  Computer  graphics  and  interactive  processing  play  a 
central  role  in  the-  activities . 

Many  commercial  organizations  are  developing  and  marketing  packages 
aimed  at  finite  element  manipulation.  Some  of  these  aim  at' enhancing  a 
specific  analysis  program.  Others  are  aimed  at  a  specific  hardware 
product,  but  address  popular  finite  element  programs.  The  third  group 
of  programs  are  both  hardware  and  analysis  program  independent.  There 
are  obvious  advantages  and  disadvantages  to  each  alternative.- 


Apart  from  the  development  of  user-oriented  ore-  and  postprocessor 
software,  research  efforts  have  been  directed  at  some  fundamental 
points  associated  with  mesh  generation  and  software  engineering.  Among 
the  subjects  addressed  are  algorithms  for  geometry  definition, 
discretization  of  irregular  shapes,  and  a  new  and  promising  subject, 
grid  optimization.  This  contribution  will  address  the  subjects  of 
geometry  definition,  mesn  generation  algorithms  available  as  part  of 
postprocessing  software,  postprocessing  and  digitizing  techniques, 


mesh  optimization,  and  the  design  o£  relevant  software.  Mo  claim  to 
completeness  is  being  made. 


In  going  beyond  analysis  to  the  associated  automatic  model 
generation  and  result  presentation,  the  process  changes  from  that  of 
"computer  analysis"  to  "computer-aided  analysis".  An  increase  in  the 
computer  costs  over  what  is  necessary  for  the  analysis  itself  is  to  be 
expected,  but  the  resulting  savings  in  labor,  and  the  faster  response, 
usually  more  than  compensate  for  the  additional  expense. 

The  process  of  model  generation  is  by  no  means  divorced  from  that 
of  computer  aided  design  and  computer  aided  manufacturing.  It  has  long 
been  recognized  that  an  integration  of  design  and  analysis  functions 
would  represent  an  ideal  working  environment.  Although  this  is 
indisputable  in  principle,  its  practical  implementation  is  not 
straightforward  by  any  means,  and  attempts  to  construct  a  formal  link 
between  a  design  data  base  and  analysis  programs  for  complex  objects 
have  been  typically  costly  and  laborious.  It  is  in  the  nature  of 
complex  systems  that  an  increase  in  size  brings  with  it  a  decrease  in 
reliability  and  flexibility.  We  may  still  see  some  success  in  such 
integration  attempts,  but  an  acceptable  solution  will  only  be  possible 
after  significant  advances  are  made  in  computer  software  design 
technology,  and  more  effective  use  is  made  of  current  hardware. 


Geometry  definition 


The  process  of  model  definition,  in  preparation  for  a  finite 
element  analysis,  comprises  several  distinct  stages.  First  a 
mathematical  definition  of  the  geometry  is  accomplished,  and  then  the 
discretization  process  may  be  carried  out.  In  creating  the  geometrical 
model,  mathematical  methods  for  the  representation  of  curves  and 
surfaces  in  three  dimensions,  typical  of  those  employed  in  computer 
aided  design  may  be  employed.  Some  of  the  methods  used  in  the 
description  of  surfaces  for  analysis  purposes  are  based  on  mathematics 
used  in  finite  element  computations.  A  typical  example  is  to  be  found 
in  isoparametric  element  formulation.  Such  methods  are  usually  not  as 
sophisticated  as  those  used  in  surface  design,  producing 
discontinuities  in  slopes  and  curvatures  at  patch  boundaries.  They  are 
nontheless  popular  due  to  the  familiarity  of  the  analysts  with  their 
formulation,  the  associated  ease  of  implementation,  as  well  as  their 
ability  to  describe  discontinuous  geometry.  Although  these  methods  are 
generally  adequate,  it  has  been  found  in  certain  applications,  such  as 
the  analysis  of  arbitrary  shells,  chat  the  results  may  be  extremely 
sensitive  to  curvatures,  thus  encouraging  the  use  of  the  more 
sophisticated  formulations  for  these  purposes. 
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Curves  in  three  dimensional  space  may  be  defined  ir.  parametric 
form(H.  For  curves  of  gceatar  complexity,  a  piece-wise  continuous 
parametric  representation  is  necessary.  Of  the  many  possibilities  for 
parametric  representation,  the  3ezier(2]  curves  may  be  mentioned.  The 
Bezier  method  of  surface  representation,  originally  developed  for  the 
French  automobile  industry,  derives  the  parametric  shape  of  a 
particular  curve  from  a  so-called  control  polygon,  connecting  n+1 
vertices.  The  resulting  parametric  curve  of  order  n  passes  through  the 
two  end  points,  and  is  tangential  to  the  polygon's  first  and  last 
sides.  The  shape  of  the  control  polygon  influences  the  shape  of  the 
curve.  The  control  exerted  by  the  polygon  is  global  in  the  sense  that 
a  change  of  the  position  of  one  of  the  points  will  produce  a  change  of 
the  shape  of  the  curve  overall.  If  relatively  few  points  are  used  in 
the  polygon,  the  control  of  the  shape  of  the  curve  is  good.  A  Bezier 
curve  is  a  surface  obtained  by  forming  the  product  of  two  3ezier 
curves . 

The  use  of  3ezisr  method  presents  some  problems  in  the  control  of 
the  shape  of  the  curves  and  surfaces,  in  addition  to  problems  of 
geometric  continuity.  Another  approach  which  avoids  these  shortcomings 
is  that  based  on  spline  functions.  Such  functions  involve  the  fitting 
of  a  piecewise  continuous  curve  through  a  number  of  points,  and 
provides  for  continuity  of  geometric  derivatives  at  segment  boundaries 
through  constraint  conditions.  It  has  local  control,  so  that  the 
effect  of  moving  a  point  on  the  shape  of  the  curve  applies  only  to  the 
immediate  neighbourhood  of  the  point.  It  is  possible  to  form  closed  or 
open  curves,  as  in  the  case  of  the  3ezier  method,  as  well  as  introduce 
corners.  Cubic  splines  are  the  ones  most  often  used. 

There  exists  a  number  of  works  which  address  the  question  of  curve 
and  surface  description,  and  ref.  (11  could  be  consulted  for  additional 
decail.  For  example,  Rawat (31,  describes  the  use  of  bivariate  splines 
in  the  definition  of  ship  hull  geometry,  ^hailhe imer  and  Mckee[4] 
survey  the  use  of  spline  functions  'in  the  description  of  ship  hull 
surfaces.  Among  the  topics  covered  are  the  construction  of  complex 
curves  from  control  polygons,  the  computation  of  surface  intersections 
and  the  solution  of  the  hidden  line  problem  associated  with  such 
3ur faces . 

Gordon(3j  develops  a  blended  spline  approach  to  define  complex 
surfaces  using  oararaecric  boundary  curves.  An  example  of  the  use  of 
cubic  splines  as  a  basis  for  a  mathematical  surface  definition, 
followad  by  a  finite  element  analysis,  is  given  by  Liu(5].  He 
describes  an  approach  by  which  the  smoothness  of  the  surface  generated 
for  arbitrary  curved  shell  analysis  nay  be  optimized  using  the  unknown 
twist  vectors  at  the  patch  corners  as  controlling  parameters. 
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Mash  generation  algor  itluns. 


Many  papers  exist,  which  deal  with  mesh  generation  algorithms,  soe 
for  example  [7,3,9,10],  an  early  survey  of  mesh  generation  algorithms 
is  found  in  [31].  In  this  section  we  deal  primarily  with  methods  for 
the  automatic  subdivision  of  complex  geometrical  configurations,  and 
the  associated  computation  of  the  position  of  intermediate  grid 
points.  The  simplest  and  most  widely  solved  problem  is  the  two 
dimensional  one,  in  which  a  complex  shape,  or  an  interconnected  set  of 
arbitrary  shapes,  is  subdivided  into  triangular  or  quadrilateral 
subdomains.  Next  in  degree  of  difficulty  is  the  problem  of  definition 
and  discret  izat  ion  of  three  dimensional  surfaces,  which  may  be 
in terconnec ted .  The  additional  problems  arise  from  the  need  to 
determine  a  method  for  the  definition  of  the  location*  of  the  surface 
in  space,  as  well  as  the  fact  that  several  surfaces  may  intersect 
along  the  same  curve.  At  such  intersections  surface  parameter 
continuity  may  be  required,  although  not  necessarily.  Finally  we  have 
the  case  of  three  dimensional  solid  continua,  where  a  volume,  bound  by 
such  surfaces,  may  be  subdivided  into  a  number  of  tetrahedron  or 
be  ick-type  el  aments . 

It  has  been  mentioned  before  that  model  creation  involves  both 
mathematical  definition  in  continuous  ‘form,  and  subsequent 
discretization  into  finite  elements  (subspaces).  The  procedure  by 
which  this  may  be  accomplished  is  not  unique,  and  may  take  several 
forms,  as  is  evident  from.  figura(Lj.  It  is  possible,  for  example  to 
generate  continuous  forms  for  both  curves  and  surfaces,  and  then 
discretize  both  independently .  Since  the  curves  form  the  surface 
boundaries,  care  must  be  taken  that  the  discretization  is  compatible. 
Alternately,  one  may  proceed  from  a  continuous  curve  directly  to  a 
discrete  curve  format,  and  use  this  information  to  form  the  discrete 
surface  form  without  necessarily  providing  a  continuous  surface 
definition.  It  is  also  possible  to  define  curves  in  dicreta  form 
directly,  say  by  digitizing  point  by  point,  and  then  proceed  to 
generate  discrete  forms  for  both  surfaces  and  solids. 

Another  approach  evolves  from  disregarding  the  natural  geometrical 
hierarchy,  and  proceeding  directly  to  th  2  definition  of  the  surfaces 
in  continuous  form.  Once  the  surfaces  have  been  defined,  their 
in ters.ee cion  is  computed,  thus  defining  their  boundary  curves,  as  well 
as  the  form  of  the  solids  bounded  by  them.  Ail  the  above  approaches 
have  their  advantages  and  disadvantages .  A  truely  general  purpose 
model  generation  program  should  allow  all  possible  oaths,  so  that  the 
user  may  choose  whichever  approach  may  best  suit  his  particular 
ocoblem . 


Turning  to  mesh  generation  algorithms,  a  method  for  subdividing  2 
two  dimensional  area,  bound  by  an  arbitrary  polygon,  into  triangular 
finite  elements  is  given  by  3ykat(ll].  The  method  subdivides  the  area 
first  into  a  number  of  convex  subregions.  This  is  followed  by  an 
automatic  "pre-conditioning"  of  the  boundaries,  by  introducing 

additional  boundary  points  to  eieminate  sudden  changes  in  grid 
spacing.  Automatic  boundary  refinement  is  applied  at  re-entrant 
corners.  Finally  all  convex  regions  are  triangulated,  producing  well 
shaped  elements.  A  two  dimensional  mesh  generator  with  many 

interesting  features  is  described  by  Gabr ieison [ 1 5 ] .  It  is  based  on 
the  definition  of  multi-connected  four  sided  regions  (i,j  grids),  in  i 
which  intermediate  points  are  computed  by  linear  interpolation.  It  i 
accepts  input  alphanumer ically  via  the  terminal,  in  the  form  of  data  I 
cards,  digitizer  output  or  directly  from  a  larger  geometric  data  base 
intended  for  computer  aided  manufacturing.  The  generated  meshes  can  be 
used  for  either  heat  transfer  or  stress  analysis. 

Gordon  and  Hall (121  address  the  problem  of  domain  mappings  in  the 
context  of  finite  element  analysis  and  present  a  heuristic  method. 
Nielson  and  Wixora(13]  discuss  different  techniques  for  surface  i 
definition  from  boundary  curves.  Extensions  to  the  method  of  Gordon 
are  presented.  A  mesh  generator  based  on  3-spline  curve 
representation [14]  uses  a  command  language  for  the  generation  of 
points,  lines,  circles,  and  various  other  curve  and' surface  types.  It  ■ 
has  the  ability  to  calculate  surface  intersections. 

An  early  example  of  mesh  generators  for  solid  analysis  is  described  j 
in  [19],  There,  the  finite  element  isoparametric  formulation  is  used  ; 
as  a  basis  for  surface  representation  and  discretization  as  well  as 
volume  discretization-  Cook  [23,21]  develops  and  uses  the  concept  of 
"natural  coordinates"  to  generate  surface  and  volume  grids.  Another 
effort  based  on  the  same  concept  is  "presented  by  Aral(22].  A  process  ; 
for  checking  solid  models  by  means  of  computing  and  displaying  the 
intersections  of  one  elements  of- the  model  with  a  series  of  planes  is 
descr ibed  in  [ 23] . 

Some  authors  present  attempts  at  developing  a  language  to  be  used 
as  a  means  for  input  description.  For  example,  the  SAIL  language [ 24] , 
developed  by  3oeing  is  a  FORTRAN-1  ike  structural  input  language,  which 
generates  input  data  for  many  finite  element  programs,  particularly 
NAS TRAN  and  SAM£S .  Silvester  [  25  ]  describes  the  MAG_NE7  language’,  which 
is  a  tool  for  the  description  of  finite  element"*  models.  It  allows 
geometrical  and  topological  changes,  unions  of  separate  models  as  well 
as  addition  and  removal  of  model  constituents.  Input  may  be  through  a 
keyboard  or  graphics  cursor.  It  is  designed  particularly  for 
repetitive  structures. 

In  most  programs  a  bandwidth  or  wavefront  optimizer  is  necessary  to 
reduce  subsequent  solution  times,  and  sometimes  to  avoid  exceeding 
program  limitations.  A  method  for  renumbering  elements  in  order  to 
achieve  improved  storage  in  preparation  for  a  frontal  solution,  and 
other  methods  not  requiring  a  stiffness  assembly,  is  given  by 
3ykat(25].  It  represents  a  modification  of  the  Cuthill-McXee 
algor ithm  ( 27 ] ,  and  requires  a  reasonable  choice  of  the  first  element. 
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Gray  at  al.(23i  describe  two  methods  for  automatic  renumbering  of 
nodes  for  water  network  analysis.  The  method  can  be  extended  to  finite 
element  programs.  It  is  aimed  at  the  minimization  of  matrix  fill 
during  the  solution  process. 

Other  special  efforts  deserve  attention.  An  example  is  an  algorithm 
by  ?reiss[291,  which  checks  topological  consistency  of  plate  and  solid 
models.  It  is  part  of  an  in-house  NCNSA?  preprocessor.  In  problems 
where  model  geometry  changes  due  to  chemical  reactions  or  physical 
phase  changes,  such  as  models  with  abelating  boundaries,  the  problem 
of  automatic  model  generation  is  further  complicated.  Weeks  and 
Cost  [33]  propose  an  algorithm  for  cracking  such  boundaries,  in 
conjunction  with  an  automatic  mesh  generator,  and  a  finite  element 
program  for  the  calculation  of  transient  thermal  or  stress  analysis  of 
such  sructures. 

Model  generation  programs  may  be  written  as  general  purpose 
packages  or  as  special  purpose  programs.  A  general  purpose  mesh 
generator  is  more  versatile,  being  abie  to  handle  a  wide  variety  of 
problems,  and  respond  well  to  new  and  ur.forseen  situations.  On  the 
other  hand,  it  will,  by  necessity,  be  more  complex  and  will  require  a 
certain  amount  of  training  before  it  can  be  used  affectively.  A  truelv 
general  purpose  program  should  impose  no  restrictions  on  problem  size, 
or  type.  If  required  to  run  interactively,  it  will  typically  be 
restricted  to  a  small  amount  of  core.  This  combination  of  requirements 
may  result  in  extensive  disk  usage  and  computer  overhead  to  perform 
the  necessary  I/O  operations.  The  other  alternative  is  that  of  a 
special  purpose  mesh  generator,  subject  to  certain  restrictions 
regarding  problem  type  and  perhaps  size.  This  option  is  particularly 
attractive  in  situations  where  a  certain  ty pe  of  structure  is 
repeatedly  analyzed  with  different  geometric  parameters.  In  such  a 
case  it  is  often  possible  to  produce  an  optimized  grid  based  on  the 
overall  characteristics  of  the  structure  in  mind,  and  thereby  relieve 
the  user  of  detailed  decisions  and  input  detail.  It  is  often  possible, 
and  indeed  most  desireable,  to  use  an  approach  combining  the  two 
possibilities  to  best  advantage.  In  this  case  a  trueiy  general  purpose 
mesh  generator  is  coupled  with  several  special  "pre-preprocessors'', 
whose  functions  would  be  to  provide  appropriate  instructions  for  the 
generation  of  specific  structural  models  of  known  configuration  using 
the  general  purpose  generator  and  a  minimum  of  input. 

An  excellent  example  of  a  special  purpose  mesh  generation  program 
is  that  developed  by  Leick  and  ?otvin(15].  A  special  purpose  analysis 
program  was  written  to  handle  tubular  joints  such  as  T,y,X  and  TX 
joints,  prevalent  in  offsnore  tubular  structures.  The  associated  mesh 
generator  produces  a  fine  grid  near  stress  concentration  areas  and  a 
course  one  away  from  them.  The  grid  characteristics  are  defined  in 
terms  of  a  few  parameters. 
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Available  Mesh  Generation  Software 


In  general  mesh  generators  may  be  developed  for  specific  hardware/ 
or  may  be  intended  to  be  hardware  independent.  Hardware  dependent 
packages  are  usually  associated  with  turnkey  systems,  including  tor 
example  refresh  type  graphics.  In  the  same  manner,  pre-  and 

pos tprocesso r s  may  be  tailored  to  a  specific  analysis  program,  or  may 
be  intended  for  use  by  a  number  of  different  programs.  There  are 

advantages  and  disadvantages  to  system  dependent  packages  and  system 
independent  ones.  Whether  the  package  is  dependent  on  hardware  or 
analysis  program  or  both,  it  is  usually  well  maintained  by  the 
originators.  Should  the  hardware  or  analysis  program  change,  the 
auxiliary  software  is  automatically  changed  and  tested.  Furthermore, 
the  ore-  and  postprocessor  programs  would  cater  well  to  the 
ideosynchr onac ies  of  the  analysis  program,  and  would  utilize  special 
hardware  features  most  effectively.  Cn  the  other  hand,  the  user's 
decendence  on  a  single  source  hampers  his  freedom  or  action.  Should  a 
oarticuiar  tyoe  of  hardware  be  dropped  by  the  manufacturer,  or  proves 
to  be  ineffective,  the  required  change  over  to  another  system  is 
usually  costly  in  both  lost  production  time  and  manpower.  Furthermore, 
if  the" auxiliary  software  is  analysis  program  dependent,  and  the  user 
has  several  in-house  programs,  the  personnel  training  problems  may 
become  exessive.  In  the  case  of  general  purpose  auxiliary  software, 
the  training  may  be  more  demanding,  but  need  only  be  done  once  tor  all 
programs.  A  maximum  flexibility  of  action  is  guaranteed.  Cn  the  other 
hand,  it  is  more  time  consuming  to  maintain  a  hardware  and  analysis 
orogram  independent  program.  Any  change  in  hardware  or  analysis 
program  character  is  tics  has  co  be  carefully  studied  as  to  _its  effect 
on  the  ora-  and  postprocessor  software.  The  group  performing  t.oe 
main  tai  nance  lias  to  have  continuous  access  to  a  l*  hardware  and 
software  in  question.  This  is  only  possible  for  extremely  Large 
organizations.  It  .tray  also  be  accomplished  through  a  well  organized 
users  group.  A  survey  of  available  ore-  and  postprocessor  software  is 
given  by  Mapolitano  at  al.[3  2] 

Turning  to  specific  software  reported  in  the'  literature,  it  is 
perhaps  appropriate  to  start  by  mentioning  those  packages  supporting 
popular  analysis  software.  *3  start  by  surveying  auxiliary  software 
specif  icaiiy*  a  Load  at  particular  analysis  program,  following  that,  we 
snail  consider  analysis  program  independent  software. 

Of  ail  commercial  programs,  perhaps  the  one  most  used  in  practice 
is  MASTRAN.  An  early  survey  of'  NA5TRAN  ore-  and  postprocessors  is 
given  in  ref  [IT]. The  paper  lists  a  number  of  model  generation  ar.d 
bandwidth  node  resequencing  programs  available  through  various 
organizations . 

As  mentioned  earlier,  it  is  impc 
generation  geometry  and  logic  for  later 
boundary  conditions.  A  pioneering  affori 
by  Coo<[i'3],  who  used  spline  represents' 
geometry  definition.  The  resulting  orepi 
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NA3TRAN  program.  A  general  purpose  pos tprocesso r  for  NASTRAN,  which 
operates  in  booh  the  batch  and  time-sharing  environments  is  described 
in  ref. [33].  A  modelling  system  designed  for  NASTRAN  is  described  in 
[341.  It  is  based  on  a  hierarchial  representation  of  gemoetry,  here 
called  "ccnstruction-in-contaxt" ,  and  has  parametric  representation  of 
curves , sur faces  and  volumes.  IGFES[35]  is  an  interactive  NASTRAN  pre- 
and  postprocessing  system,  operational  on  I3M-363  and  PDP-11 
computers.  It  is  modular  and  uses  storage  tube  graphics  as  well  as 
mechanical  plotters.  It  includes  both  two  and  three  dimensional  shell 
mesh  generators  and  model  editors.  Amongst  mesn  generation  options  are 
Laplacian,  linear  interpolation  between  opposite  sides  (here  called 
direct  ray) ,  Coons  blending  and  isoparametric  mapping.  It  has 
contouring  capabilities. 

Several  auxiliary  software  packages  have  been  designed  with  the 
popular  SAP  program  in  mind.  One  of  the  preprocessors  is  due  to 
Raid j ian  [  35 ] .  It  is  interactive,  efficient  but  has  somewhat  limited 
capabilities.  FE?SA?[37]  is  an  optimized  version  of  the  SAP  program 
for  CDC  computers.  It  has  specialized  two  dimensional  and  solids 
preprocessors.  These  preprocessors  may  also  be  used  to  generate  loads 
and  boundary  conditions.  German  iso her  Lloyd(33]  created  a  system  of 
programs  built  around  the  SAP- IV  program.  Among  the  available  programs 
is  an  automatic  mesh  generator  called  GLGEN. 


Two  other  programs,  particularly  popular  in  Europe  are  SESAM-S9[39] 
and  ASXA(43],  Haugecud [ 4 1 ]  describes  a  solids  preprocessor  written  for 
SESAM-S9.  It  is  based  on  the  creation  of  a  topological  arrangement  of 
three  dimensional  grids,  which  can  then  be  distorted  to  fit  the 
physical  geometry  of  the  object  being  analyzed.  Automatic  node 
numbering,  definition  of  physical  properties  and  specifications  of 
loads  are  possible.  A  further  development  of  the  3ESAM-69  program  is  a 
modular  system(42],  which  includes  a  finite  element  ore-  and 
postprocessor  utility  package,  DA5A[43).  DA3A  has  approximately  3033 
statements  and  233  subroutines.  It  required  4  to  5  man  years  to 
develop.  It  allows  the  user  to  implement  special  and  general  purpose 
data  generators,  using  its  subroutines  as  building  blocks.  It  is  based 
on  the  concept  of  the  "mesh" ,  which  is  defined  by  a  system  of 
orthogonal  surfaces,  which  may  assume  prismatic,  cylinder ioai  or 
spherical  shapes.  The  generated  medal  may  be  modified  later  to 
introduce  local  detail.  The  system  also  allows  for  automatic  mesn 
refinement,  as  well  as  load  and  boundary  condition  generation. 


Uerkvitz [ 44 j  and  Gr  ieger [  4 5 , 4 6 ]  describe  a  geometrically 
hierarchial  mesh  generation  program,  called  INGA,  aimed  at  the  finite 
element  program  A3XA.  It  uses  finite  element  interpolation  functions 
as  a  basis  for  its  mesh  generator,  and  is  interactive  in  nature,  using 
primarily  a  refresh  graphics  stand-alone  system. 


Moving  on  to  program  independent  auxiliary  software,  an  earl’ 
example  is  given  by  the  work  of  Sousquet  and  Ta tes [73] .  The  previewing 
and  model  editing  software  was  intended  as  an  analysis  orocra: 
independent  package,  and  operated  on  storage  tube  terminals,  thus 
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domain.  The  UNISTRUCJ  47]  package  uses  the  so  called  "drag  method"  in 
order  to  develop  three  dimensional  node  and  element  arrangements  from 
two  dimensional  ones,  and  two  dimensional  ones  from  one  dimensional 
strings.  The  FASTDRAW  prog  ram (43]  has  been  developed  by 
McDonnell-Qouglas  and  serves  as  a  preprocessor  for  STRUDL,  NASTRAN  and 
AN SYS .  Tektronix  has  developed  a  software  package  for  its  line  of 
graphics  hardware  devices,  called  F EM .  It  is  designed  for  a  stand 
alone  intelligent  device  such  as  the  Tektronix  FSM-131  system.  It  has 
extensive  interactive  editting  capab il ities ,  and  a  new  mesh  generation 
module  has  been  announced. 

Several  packages  of  ?re-  and  postprocessing  software  have  been 
developed  in  a  University  environment  as  a  byproduct  of  research 
efforts,  or  out  of  necessity  to  facilitate  access  to  available  finite 
element  software.  The  GIFTS  system [33]  is  one  example.  It  has  two  and 
three  dimensional  model,  load  and  boundary  condition  generation 
capabilities.  It  is  designed  to  run  in  minimum  core  space.  It's 
postprocessor  handles  all  but  solid  models  at  present.  It  is  supported 
by  a  users  group,  and  is  distributed  out  of  the  University  of  Arizona. 
Interfaces  exist  with  the  SAP-IV  and  ANSY5  programs,  and  others  are 
under  development.  Apart  from  ore-  and  postprocessing  modules,  GIFTS 
is  also  capable  of  static  and  dynamic  analysis  on  both  a  number  of 
mainframes  and  minicomputers.  A  finite  element  preprocessor, 
FEMGEN [51,53] ,  has  been  under  development  at  the  University  of  Lund  in 
Sweden-  since  1974.  It  is  intended  to  support  analysis  programs  such  as 
ASKA,  NASTRAN ,  STRUDL,  SAP-IV  and  ACINA.  FEMGEN  uses  storage  tube 
terminals,  drum  plotters  and  line  printers.  Straight  lines  and 
circular  arcs  may  be  defined  in  terms  of  geometric  points,  or  as  the 
intersections  of  standard  mathematical  surfaces,  such  as  planes, 
cylinders  and  cones.  The  program  runs  in  25X  words  on  a  CNIVAC-1133, 
and  has  been  converted  to  GDC,  IBM  and  PDP-11  computers. 


Many  other  efforts  are  underway.  For  example,  the  SURERNnT 
program(52]  is  being  developed  in  the  German  Federal  Republic  by  3raun 
3overi  {3BC).  The  program  supports  line,  surface  and  solid  elements, 
and  follows  a  geometric  hierarchy.  Interfaces  are  being  written  for 
analysis  programs  such  as  ASX A ,  NASTRAN,  MARC,  TCPAS  and  ANSY3 . 
SUPERNET  may  be  run  in  batch,  or  in  an  interactive  mode.  In  the  latter 
case  a  Computer  Vision  CAD  system  is  used.  A  ore-  and  pos tor ocesso r 
package  called  FEMALE (531  is  intended  as  a  general  purpose  package. 
However  it  seems  that  its  mesh  generation  capabilities  are  still  at  an 
early  stage.  From  the  current  description  it  appears  to  be  more  of  a 
model  pre-viewer  rather  than  a  preprocessor.  An  interactive  three 
dimensional  shell  mesh  generator  is  described  in  ref (54]  .  It  uses 
refresh  scope  graphics,  light  pen  and  keyboard.  Initial  input, 
defining  the  outlines  of  surface  patches  is  introduced  via  cards.  The 
user  may  view  the  outlines  and  perform  certain  operations.  Afterwards, 
the  mesn  generator  may  be  invoked,  and  results  displayed  and  editted. 
A  bandwidth  optimizer  is  included,  and  output  may  be  in  the  form  of 
punched  cards,  permanent  files  or  magnetic  tape  output.  Crawford[55J 
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is  in  the  process  of  developing  a  pre-  and  postprocessor  package, 
which  uses  a  geometrical  hierarchy. 

Some  programs  have  an  integrated  pre-  and  postprocessing 
capability,  which  may  often  be  important  in  special  purpose  codes,  in 
which  the  processes  of  model  generation  and  of  computation  are 
interleaved.  An  interresting  example  is  the  TOTAL  program  reported  cy 
3eaubien [ 57J .  It  is  a  two  dimensional  analysis  program  designed  to 
simulace  the  failure  mechanisms  in  layered  orthotropic  composites.  It 
has  a  built  in  pre-  and  postprocessor  which  helps  keep  track  of  the 
progress  of  cracks. 

Qualitative  evaluation  of  interactive  versus  batch  computations  are 
hard  to  find  in  the  literature.  Measurements  are  bound  to  be  system 
dependent,  and  often  represent  personal  opinions  and  experiences.  It 
is  nevertheless  interesting  to  mention  figures  extracted  from  a 
promotional  pub!ication[34]  .  It  is  cla-imed  there  that  an  interactive 
solution  costs  only  45%  of  the  cost  of  a  similar  computation  performed 
in  the  batch  mode.  Gut  of  133  dollars  spent  in  a  batch  computation,  73 
are  scene  on  modal  generation,  13  on  computer  processing  and  23  on 
result  interpretation.  In  the  interactive  solution,  45  dollars  are 
spent  on  the  computation,  23  of  which  are  spent  on  model  generation, 
15  on  computer  processing  and  13  on  result  interpretation.  It  is  not 
clear  whether  these  figures  are  typical  of  the  program  in  question  or 
are  generally  valid. 


Postprocessing  and  display  techniques. 


Although  postprocessing  is  often  directly  related  to  the  mesh 
generation  process,  special  display  techniques  are  often  employed  at 
this  stage,  so  that  it  is  appropriate  to  report  on  postprocessor 
systems  and  such  techniques  under  the  same  heading. 

Contouring  of  stress  and  deflection  results  is  a  useful  and  popular 
postprocessing  operation.  Re  f . [ 5  3 ] ,  for  example,  gives  a  "CRT RAN 
program  which  admits  several  subregions,  each  divided  into  triangular 
or  quadrilateral  finite  elements.  The  problem  of  contouring  within 
higher  order  element  boundaries  is  mathematically  intriguing.  An 
initial  effort  has  been  conducted  by  Akin  and  Grav(59i,  in  which  they 
lay  out  the  fundamental  relationships  and  provide  some  numerical 
examples.  A  follow  up  paper [53]  utilizes  a  predictor-corrector  method 
to  trace  suc.n  contours  accurately.  Similarly,  Meek  and  Beer  [51] 
describe  a  method  for  tracing  deflections  across  higher  order 
elements.  It  is  based  on  a  process  involving  averaging  .nodal  values, 
followed  by  element  subdivision  and  local  linearization.  Results  are 
given  for  triangular  and  quadrilateral  elements. 
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Akin  and  3toddart[S2]  peasant  an  algorithm  for  plotting 
intersections  of  isooar amatr ic  solid  elements  with  an  arbitrary  plane, 
and  tne  subsequent  plots  of  strass  or  displacement  contours  on  the 
resulting  mesh.  The  plots  are  generated  point  by  point  in  a  random 
order,  and  are  suitable  for  graphics  terminals  or  electrostatic 
plotters.  Frey  ec  al .  [53]  describe  a  method  for  calculating  the 

intersection  of  a  plane  with  second  ordsr  isoparametric  solid  brick 
elements . 


A  paper  describing  an  algorithm  for  clipping  and  capping  solid 
polyhedron! 54 ]  ic  directly  applicable  to  the  display  of  three 
dimensional  models.  It  produces  an  input  file  for  the  well  distributed 
hidden  line  program,  SCVIS-3VU [ 55 ] .  A  version  of  the  program  foe  15 
bit  minicomputers,  called  MC  VI H  -  A  K I Z  C  A ,  is  also  available  [35]  .  A  two 
and  three  dimensional  half-tone  and  color  program  foe  the  output  of 
finite  element  results  is  described  by  Murai  and  Tataishi [55] .  It  has 
some  mesh  generation  capabilities,  is  capable  of  solving  the  hidden 
line  and  surface  problems  and  contouring.  Furthermore  it  uses 
inexpensive  hardware. 


Mesh  Optimization 


Among  the  most  interesting  areas  of  current 
preprocessing  is  that  of  mesh  optimization.  Here  the 
devising  the  optimum  mesh  configuration  for  a  particular 
and  an  approximate  number  of  degrees  of  freedom  in  order 
optimum  solution  for  the  problem  at  hand.  It  is  obvious 
definition  of  what  constitutes  an  optimum  solution  is 
will  be  found,  for  example  that  high  accuracy  in  a  peak 
achieved  by  an  optimized  grid  may  be  associated  with  poor  deflection 
results.  It  is  also  obvious  that  a  grid  may  be  optimum  for  one 
specific  loading  case,  but  far  from  it  for  another. 
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Several  authors  have  been  active  in  this  area,  as  is  evident  from 
r ef  .  [  57 , 53  ,  5 9 , 73 , 71 ] .  A  paper  by  Turks  and  McNeicef72]  introduces  the 
node  coordinates  of  the  mesh  as  independent  unknowns  in  the  finite 
element  variational  formulation.  The  solution,  which  is  based  on  a 
minimization  of  the  potential  energy  of  the  system,  provides  both 
displacement  values  as  well  as  new  node  coordinates.  Certain 
constraints  have  to  be  imposed  on  the  node  relocation  to  ensure 
meaningful  results.  Carroll  [73]  extends  grid  optimization  to  the  area 
of  vibrational  mode  computation.  He  shows  that  a  different  optimum 
grid  is  associated  with  each  mode. 
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Turke  and  hcNeice[74]  and  Meiosh  and  Xiiiian[7!5]  show  the 
superiority  of  an  optimized  mesh  over  one  obtained  by  uniform 
refinement.  The  introduction  of  additional  freedoms  and  elements  was 
allowed  in  [75],  based  on  local  solutions  to  study  the  effect  of 
element  subdivision.  Denayer[76]  generalized  some  algorithms  given  in 
[3]  for  the  generation  of  topologically  uniform  grids.  The  mapping 
process  from  a  parent  grid  to  the  area  under  consideration  is  chosen 
so  as  to  achieve 'cer tain  optimal  conditions,  which  translate  into  a 
variational  formulation  lending  itself  to  a  finite  element  solution. 
Carey(77]  presents  a  method  for  selective  refinement  of  a  finite 
element  mesh  near  an  area  of  interest.  Lagrangian  constraints  are  used 
to  maintain  continuity  with  the  rest  of  the  model.  Results  show  good 
convergence  in  the  solution  of  the  Poisson,  equation  with  a  singularity 
in  two  dimensions. 


Interactive  Software  Design 


Associated  with  interactive  software,  special  problems  in  design 
methodology,  software  reliability,  portability  and  testing  arise.  Some 
of  the  surveyed  work  is  being  directed  to  this  specific  area. 
Requirements  for  developing  general  purpose  ore-  and  postprocessing 
systems  are  formulated  by  Tischier  and  3ernier[79],  and  discussed  by 
Berness  and  Tocher.  In  [33]  it  is  claimed  that  73%  to  93%  of  the  cost 
of  a  finite  element  analysis  is  typically  attributed  to  manpower 
costs,  the  rest  being  computer  charges.  The  concept  of  a  modular 
program  library,  in  which  the  ore-  and  postprocessors  interact  with 
the  analysis  programs  via  a  standard  data  base  is  well  represented  by 
the  GIFTS  packags[31].  The  design  of  interactive  ore-  and 
postprocessor  systems  is  described  in  [32],  which  discusses  procedures 
for  program  design,  testing,  overlaying,  data  handling  and  command 
language  design.  The  problem  of  program  portability  is  more  complex 
for  graphics  software  than  it  is  for  standard  batch  analysis  programs. 
Foley[33]  discusses  the  design  of  device  independent  basic  software. 


Conclus ion 


The  paper  reviews  the  principal  components  of  finite  element  pre- 
and  postprocessing  (so-called  "auxiliary  finite  element  software").  It 
provides  an  assesment  of  the  current  state  of  the  art  and  a  literature 
survey . 
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